Tunneling spectra for quasi-one-dimensional organic superconductors 
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The tunnehng conductance spectra of a normal-metal/insulator/quasi-one-dimensional supercon- 
ductor is calculated by using the Blonder- Tinkham-Klapwijk formulation. The pairing symmetry of 
the superconductor is assumed to be p, d, and /-wave. It is found that there is a well defined zero 
energy peak in electron tunneling along the direction parallel to the chains or normal to these when 
the transmitted quasiparticles feel different sign of the pair potential. The actual line shape of the 
spectra is sensitive to the nodes of the pair potential on the Fermi surface. 
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I. INTRODUCTION 



Almost 20 years after the discovery of the organic su- 
perconductors 1^] , the problem of determination of their 
pairing state has not yet found a definite solution. The 
critical magnetic field Hc2 exceeds the Pauli paramag- 
netic limit and indicates that the pairing symmetry is 
triplet Ql. The Knight shift does not change between 
the normal and superconducting states and it is a sig- 
nature for triplet pairing state ||3|. However the absence 
of Hebel-Slichter peak and the power- low decay of 1/Ti 
below Tc 1^ is an indication of the presence of nodes of 
the pair potential on the Fermi surface. 

The scattering theory can be used to distinguish the 
symmetry of the pair potential In rf-wave super- 

conductors the pair potential changes sign under a 90°- 
rotation. So under appropriate orientation of the a-axis 
within the ah plane of d-wave superconductor the trans- 
mitted quasiparticles feel different sign of the pair poten- 
tial. This results in the formation of bound states within 
the energy gap, which are detected as zero energy peaks 
(ZEP) in the conductance spectra ||7|"[l0|. 

The scattering theory has also been used for the deter- 
mination of the pairing symmetry in (TMTSF)2X by Sen- 
gupta et al. However their calculation is restricted 
to the presence or absence of ZEP at the surface, which 
can not distinguish the p from the /-wave case. In a 
more realistic calculation Tanuma et al. [l^ used the ex- 
tended Hubbard model on a quasi-one-dimensional lat- 
tice at quarter-filling, to study the quasiparticle states 
near the surface of a quasi-one dimensional organic su- 
perconductor, where the pairing symmetry can actually 
be distinguished from the overall line shape of the surface 
density of states (SDOS) and the presence or absence of 
ZEP. 

In this paper we extend the BTK formula to 
calculate the tunneling conductance in a normal- 
metal/insulator/quasi-one-dimensional organic super- 
conductor, where the structure of the Fermi surface is 
taken into account. The pairing symmetry of the su- 
perconductor is assumed to be triplet p, singled d, d^y, 
and triplet /-wave. In particular, we find that the ZEP 



appears in the electron tunneling along the a or 6 axis, 
when the transmitted quasiparticles experience different 
sign of the pair potential. Also the line shape of the 
spectra is sensitive to the presence or absence of nodes 
of the pair potential on the Fermi surface. The present 
calculation is useful for (TMTSF)2PF6, although this 
salt needs pressure to show superconductivity, but is 
not very suitable for (TMTSF)2C104, which is the only 
material exhibiting superconductivity in the Bechgaard 
salts at ambient pressure, but its unit cell is doubled 
due to anion ordering Also tunneling spectroscopy 
on the quasi two-dimensional organic superconductor k- 
(BEDT-TTF)2Cu(NCS)2 has been performed by using 
scanning tunneling microscopy [Q. In this compound 
the conducting Cu(NCS)2 layers run parallel to the be 
plane, along the a-axis, and the in-plane tunneling data 
strongly suggest rf-wave pairing state. These features can 
be used to distinguish the pairing symmetry in quasi-one- 
dimensional organic superconductor. 



II. THE MODEL FOR THE NS INTERFACE 

The motion of quasiparticles in inhomogeneous su- 
perconductors is described by the four component Bo- 
goliubov de Gennes (BdG) equations. The BdG equa- 
tions are decoupled into two sets of (two component) 
equations, one for the spin up electron, spin down hole 
quasiparticle wave functions (u|(r),ti|(r)) and the other 
for (u|(r), W|(r)). The BdG equation for spin index 
s{s) =T (i) or s{s) =i (t), read § 



{Heir) - siiBH)us{r) + J dr' A,s{s,x)vj{r') 
Jdr'Al^{s,x)u,{r') - {H;{r) - sy.BH)v^{r) 



EsUs{r) 
(1) 



The single-particle Hamiltonian is given by He{r) — 
—h^ Vr /2?7T-e + ^(^) " Ep, H IS the external magnetic 
field, Es is the energy measured from the Fermi energy 
Ep. For a given spin projection s the magnetic field shifts 
the energy by ^s^isEt- A.ss{s,x) is the pair potential, 
after a transformation from the position coordinates r, r' 
to the center of mass coordinate x = {r + r')/2 and the 
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relative vector s = r ~ r' . After Fourier transformation 
the pair potential depends on the relative wave vector 
k and x. In the weak coupling limit k is fixed on the 
Fermi surface (|fe| = kp), and only its direction is vari- 
able. After applying the quasi-classical approximation, 



I.e., 



Us{r) 
v-{r) 



Us[r) 



(2) 



'(2izm] -I- z^j^j) + (2izm] -I- z2j^j)n+n_ 

4 + "4x1"+''-'^-'^+ 



(7) 



where zq = ^T[i] ~ ^T^- ^"^^ BCS coherence factors 
are given by 



= [1 + Vi?2-|A±|Vi?]/2, 



(8) 



so that the fast oscillating part, of the wave function is 
divided out, the BdG equations are reduced to the An- 
dreev equations M 



EsUs{r) = -ivpk ■ \7Us{r) + A^sik, r)w-(r) 
EsVj{r) = ivpk ■ V%(^) + ^L{l«,r)usir) 



(3) 



where the quantities Us{r) and tJ-(r) are electron-like and 
hole-like quasiparticles with spin index s, and s respec- 
tively, and is the Fermi velocity. We consider the 
normal- metal/insulator/superconductor junction shown 
in Fig. 1^. The electron momentum parallel to the inter- 
face fc|| is conserved. For the interface that is normal to 
the a-axis, the insulator is modelled by a delta function, 
located at a; = 0, of the form V6{x). The temperature is 
fixed to K. We take the pair potential as a step func- 
tion, i.e., Ass(fc||,r) — Q{x)Ass{k\\). For the geometry 
shown in Fig. |^, Eqs. ^ take the form 



EsVs(x) = iVFkFx-£vs{x) + Af^(fc||)Us(x) 



(4) 



When a beam of electrons is incident from the normal 
metal to the insulator, with momentum fc, the general 
solution of Eqs. (0), is the two-component wave function 
— (u|[|],wm]) which for x < is written as 



"TUl 



^ — ixk f. 



(5) 



where , b-^n] , are the amplitudes for Andreev and nor- 
mal refiection for spin-up(-down) quasiparticles, kp^ is 

the x-component of the Fermi wave vector kp 



_ 2rnEF_ 



The wave vector of quasiparticles in the normal-metal, 
and the wave vector of the electron-like and hole-like 
quasiparticles in the superconductor are set equal. Since 
the translational symmetry holds in the y-axis direction, 
the momentum parallel to the interface is conserved. Us- 
ing the matching conditions of the wave function a,t x = 
0, */(0) = *//(0) and *^^(0)-^';(0) = i2mV/h'^)'^ i{0), 
the Andreev and normal reflection amplitudes a|[x]i^T[i] 
for the spin-up(-down) quasiparticles are obtained as 



4n4 



^T[i] 



4 + 



'TUl 



(6) 



[1 - v/i?2-|A±|Vi?]/2, 



(9) 



and n± = v±/u±. The 
the energy gap is given by 
(A_), is the pair potential 
ted electron-like (hole-like) 
According to the BTK 
the junction, am]{Es, k\\), 
cles, is expressed in terms 

^TUl'^TU] as [| 



internal phase coming from 
(/.± = [A±/|A±|], where A+ 
experienced by the transmit- 
quasiparticlc. 

formula the conductance of 
for up (down) spin quasiparti- 
of the probability amplitudes 



(10) 



The tunneling conductance, normalized by that in the 
normal state is given by 



aiE)=a^{E^)+aiiE^), 



where 



R 



dfc||[crAr^(fc||) -f crjv^(fc||)]. 



(11) 



(12) 



(13) 



(14) 



T[i] 



III. 



PAIRING STATES AND FERMI SURFACE 
LINE SHAPE 



For the spin triplet pairing state the Cooper pairs have 
spin 1 degree of freedom. The gap function is a 2 x 2 sym- 
metric matrix which in the spin space can be written as 



A(A;) = iay{d ■ ct). 



(15) 



where a denotes the Pauli matrices and d is a vector 
which defines the axis along which the Cooper pairs have 
zero spin projection. In the following we will take d || a, 

0, 



i.e., parallel to the chains. In that case A|| 
while A|x = A^t = A(A:||). 



Mi 
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The Fermi surface (FS) consists of two branches and 
is open in the ky direction, as shown in Fig. ^. 
We consider the following pairing cases 

a) In case of pa^-wave, p^-wave superconductor 

Ap^ = Ao sin(2fc^a), Apj^ = Aosin(fcya). (16) 

Ap^ changes its sign along the FS in the direction 
while the Ap^ changes its sign along the FS in the ky 
direction as seen in Figs. ||(a) and ||(d). 

b) In case of c?a;-wave, dy-wave superconductor 

Arf^ = Aq cos{2ka:a),Ady = Aq cos{kya) (17) 

Ad^ changes its sign along the FS in the k^ direction 
while the A^^^ changes its sign along the FS in the ky 
direction as seen in Figs. |^(b) and||(e). 

c) In case of /a;-wave, fy-wave superconductor 

A/^ = Aosin(4fc^a), A/^ = Aosin(2fcj,a) (18) 

Af_^ changes its sign along the FS in the k^ direction 
while the Af^ changes its sign along the FS in the ky 
direction as seen in Figs. ||(c) and ||(f). 

d) In case of da;y-wave, dj.2_j,2-wave superconductor 

Ad^^ = Aosin(2fc:ra) x sin(fcj^a), (19) 
^rf,2_„2 = Ao(cos(2A;j;a) - cos(fcj^a)), (20) 

Ad^y, Ad 2_ 2 change sign along the FS as seen in Figs. 
1(a) and 1(b). 



IV. TUNNELING CONDUCTANCE 
CHARACTERISTICS 

In Figs. we plot the tunneling conductance (t{E) 
as a function of E/Aq for various values of zq, for differ- 
ent orientations normal to the a or b axis. The pairing 
symmetry of the superconductor is Pa;-wave, Py-wave, in 
Fig. 3, c?a;-wave, dy-wave, in Fig. ||, f^-wave, fy-wave, in 
Fig. 6[ dxy-^ave, da;2_y2-wave, in Fig. |^. 

The conductance peak is formed in the electron tunnel- 
ing along the a or 6 axis, when the transmitted quasipar- 
ticles experience different sign of the pair potential. Also 
the line shape of the spectra is sensitive to the presence 
or absence of nodes of the pair potential on the Fermi 
surface. 

For the p^-'wave case, when the a-axis of the crystal 
is at right angle to the interface, the horizontal line in 
Fig. 1^ (a), representing the scattering process, connects 
points of the FS with opposite sign of the pair potential 
for —n/a < ky < -K/a. As a result a peak exists in the 
conductance spectra, at = as seen in Fig. ^ (a) for 
Zq — 2.5. The height of the ZEP is proportional to the 
range of ky for which sign change occurs and it is ex- 
pected to have its maximum value when the a-axis of the 



crystal is normal to the interface, since for this orienta- 
tion the transmitted quasiparticles feel a different sign of 
the pair potential for all — vr/a < ky < n/a. On the other 
hand, when the a axis is along the interface, then the ver- 
tical line in Fig. || (a) connects points of the FS with the 
same sign of the pair potential for — 7r/2a < k^ < 7r/2a, 
i.e., there is no kx for which, the transmitted quasipar- 
ticles feel the sign change of the pair potential, and no 
ZEP is formed as seen in Fig. I (b). For the Px-wave 
case the nodes of the pair potential do not intersect the 
FS and the line shape of the spectra is U-like as in the 
case of the s-wave superconductor. 

The situation is opposite in the Py-wave case, where 
the scattering process for the surface orientation normal 
to the a (6)-axis, described by the horizontal (vertical) 
line in Fig. || (d) connects points of the FS with the 
same (opposite) sign. As a consequence for the surface 
orientation normal to the a-axis, there is no ky for which 
the transmitted quasiparticles feel the sign change of the 
order parameter, and no ZEP is formed, as seen in Fig. 
^ (c), while for the interface along the a-axis the trans- 
mitted quasiparticles feel the different sign of the pair 
potential, for all — 7r/2a < kx < 7r/2a, and a ZEP exists 
as seen m Fig. I (d). For the Pj,-wave case the nodes of 
the pair potential intersect the FS and the spectra has a 
V-shaped form as in the case of d-wave superconductors. 

For the d-wave case, and for tunneling along the a or 
6-axis, the scattering process, described with the horizon- 
tal or vertical line in Figs. |^(b) and ^(e) connects points 
of the FS with the same sign in all cases. This means 
that the pair potential does not change sign and no ZEP 
is formed as seen in Fig. ||. Also due to the presence of 
nodes of the pair potential along the FS, the line shape 
of the spectra is V-like. 

For the /-wave case, although the pair potential has a 
different structure, the lines in Figs. ||(c) and ||(f) con- 
nect points of the FS with the same sign-change of the 
pair potential as in the p-wave case. As a result the tun- 
neling shows ZEP as in the p-wave case, due to the sign 
change of the transmitted quasiparticles. However, in 
all cases in the /-wave case the pair potential intersects 
the FS and nodes are formed. As a consequence the line 
shape of the spectra is V-like, as seen in Fig. || unlike 
to px-'wave case, where the pair potential is nodeless and 
the tunneling spectra has a U-shaped form. The conclu- 
sion is that the line shape of the spectra can be used to 
distinguish the Px from the /-wave pairing state. Our 



results are comparable to that of Tanuma et al. |12| al- 
though their calculation was based on a different model, 
i.e., the extended Hubbard model on a quasi-one dimen- 
sional lattice at quarter-filling. 

For the dxy-'^ave (d^2_y2-wave) case, and for tunnel- 
ing along the a or 5-axis, the scattering process, described 
with the horizontal or vertical line in Figs. |^(a) and ^(b) 
connects points of the FS with different (the same) sign. 
This means that the pair potential felt by the transmitted 
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quasiparticles changes (conserves) sign during the scat- 
tering process and a ZEP (no ZEP) is formed at E = 
for the dxy-wave ((i^2„j,2-wave) case as seen in Fig. |^. 
Also due to the presence of nodes of the pair potential 
along the FS, the line shape of the spectra is V-like. In 
addition the tunneling spectra for the dxy-yfave, for tun- 
neling direction along the a or 6 axis, is equivalent to the 
tunneling spectra for (i^2_y2-wave with surface orienta- 
tion tilted by 7r/4, i.e., along the nodal direction. 

In the metallic limit {zq = 0), (j{E) has the same form 
for each pairing potential independently from the orien- 
tation. 

In table | we summarize the results concerning the 
presence or absence of ZEP, and also the overall line 
shape of the tunneling spectra for various pairing po- 
tentials. 



V. BOUND STATE ENERGIES 



These features are explained if we calculate the en- 
ergy of the midgap state, which is given for large zq by 
the value in which the denominator of Eqs. (^) and (0) 
vanishes. The equation giving the energy peak level is 
written as lllOl 



1.0. 



(21) 



In the Px-wave case, for surface orientation normal to a- 
axis, this equation has the solution E = 0, for —ir/a < 
ky < Tr/a, since n+n-\E=o — —1, and also the transmit- 
ted quasiparticles feel a different sign of the pair poten- 



tial, i.e.. 



\e=o = — 1- When a midgap state exists 



the tunneling conductance as{E,ky) is equal to 2 and 
the peek in a{E) seen in Fig. ^(a), is due to the nor- 
mal state conductance i?Ar in Eq. which depends 
inversely on the Zq for large zq. For surface orienta- 
tion normal to &-axis the range of kx for which Eq. |2l| 
has solutions collapses to zero, and no bound states are 
formed. Then a{E) goes to zero as l/zg and there is no 
conductance peak. In the Py-wave case for surface orien- 
tation normal to a-axis Eq. 21 has no solutions because 
0-0+|s=o = Ij while for surface orientation normal to 
6-axis a bound state is formed for all kx in the interval 
— 7r/2a < kx < 7r/2a. For the d-wave case and for tun- 
neling along the a or 6 axis of the crystal, the condition 
(I>-4>'^\e=o = — 1 is not satisfied and no bound states are 
formed. For the /-wave case that condition is the same 
as in the p-wave case, and thus the ZEP is formed as in 
the p-wave pairing state. For the dxy-'wave (dx^-y^-'wave) 
case that condition (j)-(f>'^\E=o = — 1 is (is not satisfyied 
for tunneling along the a and b axis and hence the ZEP 
(no ZEP) is formed for both timncling directions. 



VI. MAGNETIC FIELD EFFECT 

In this section we describe the effect of the external 
magnetic field H in the direction parallel to the chains, 
in the spectra. The effect of the magnetic field depends 
on the spin, of the quasiparticles. The tunneling conduc- 
tance is given by 



a{E) ^ a^{E - hbH) + ai{E + hbH). 



(22) 



In Fig. 1^ the tunneling conductance <7{E) is plotted for a 
fixed magnetic field (IbH/Aq = 0.2, and barrier strength 
zq — 2.5. The pairing symmetry of the superconductor 
is p, d, /, dxy respectively. 

The magnetic field splits symmetrically the tunneling 
spectrum which is a linear superposition of the spectra for 
spin up(down) quasiparticles. The amplitude of the split- 
ting depends linearly on the magnetic field H. For the 
case of /isiJ/Ao — 0.2, seen in Fig. || the spin up(down) 
part of the spectra partially overlaps while for larger val- 
ues of the magnetic field the spin up and down branches 
are well separated. In the latter case the right (left) 
branch of the spectra originates from spin up(down) 
quasiparticle spectra aq^ {E — ^BH){aq^ {E + ^bH)). 

The condition for the formation of bound states is 
slightly modified under the presence of a magnetic field to 
E — (J-bH = 0, for the spin- up region, and E + (IbEI = 0, 
for the spin-down, from the corresponding _E = in the 
absence of any field. So the multiplication of the bound 
states and the presence of magnetic field results into the 
appearance of double peak in the conductance spectra. 



VII. CONCLUSIONS 

We calculated the tunneling conductance in normal- 
metal/insulator/quasi-one-dimensional superconductor, 
using the BTK formalism. We showed that the ZEP 
appears in the electron tunneling along the a or b axis, 
when the transmitted quasiparticles experience a differ- 
ent sign of the pair potential. Also the line shape of the 
spectra is sensitive to the presence or absence of nodes 
of the pair potential on the Fermi surface, and results 
to a U-shaped structure for the p^^-wave case and to a 
V-shaped one for the d and /-wave cases. 

The ZEP are due to the formation of bound states 
within the gap. The calculation of the conductance 
WsiE, ky), for which bound state occurs, shows an en- 
hancement at the bound state energy. The effect of the 
magnetic field in the direction parallel to the chains is to 
split linearly the ZEP. 

Throughout this paper the spatial variation of the or- 
der parameter near the surface, which depends on the 
boundary orientation, is ignored for simplicity. However, 
since the features presented here are intrinsic which are 
generated by the existence of surface bound states, the 
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essential results do not change qualitatively. Also we as- 
sumed perfectly flat interfaces in the clean limit, so any 
impurity scattering and the effect of the surface rough- 
ness are ignored. 
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TABLE L Presence (Yes) or absence (No) of a zero-bias 
conductance peak in electron tunneling along the a and h axes 
for different symmetries of the pair potential. Also the overall 
line shape of the tunneling spectra is shown 



Symmetry a-axis ZEP 6-axis ZEP line shape 

pi Yes No U 

Py No Yes V 

No No V 

dy No No V 

U Yes No V 

fy No Yes V 

dxy Yes Yes V 

d^2_„2 No No V 
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(a) 

Normal Metal 



Quasi ID Superconductor 



electron 




hole-like 
quasipaiticle 



(b) 

hole-like Quasi ID Superconductor 
quasiparticle 

electron-like 
quasiparticle 




Normal Metal 

FIG. 1. The geometry 

of the normal-metal/insulator/quasi-one-superconductor in- 
terface. The arrows illustrate the transmission and reflection 
processes at the interface, (a) The insulator (vertical line) is 
normal to the a-axis. (b) The insulator (horizontal line) is 
normal to the 6-axis. 



b) 




0.5 -0.5 




0.5 -0.5 



d) 



FIG. 2. The shape of the Fermi surface (solid line) for 
the pair potentials (a) pa;-wave, (b) d^j-wave, (c) fx-wa,ve, (d) 
Py-wave, (e) d^-wave, (f) /^-wave. Inside the shaded (white) 
region the pair potential is negative(positive). kx, ky axe in 
units of 7r/a, a is the crystal lattice spacing. The horizontal 
(vertical ) arrow indicates the change in the momentum in 
the electron tunneling along the a (b) axis. 
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FIG. 3. The shape of the Fermi surface (solid line) for 
the pair potentials (a) d^y-wave , (b) d^2_y2-wave. Inside the 
shaded(white) region the pair potential is negative(positive). 
k:c,ky are in units of n/a, where a is the crystal lattice spac- 
ing. The horizontal (vertical ) arrow indicates the change in 
the momentum in the electron tunneling along the a (b) axis. 
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E/A„ 

FIG. 4. Normalized tunneling conductance cr{E) as a 
function of E/Ao for = (solid line), zo = 2.5 (dotted 
line), for different orientations; (a) and (c) for surface orien- 
tation normal to a-axis, (b) and (d) for surface orientation 
normal to b-ajcis. The pairing symmetry of the superconduc- 
tor is px in (a) and (b) and py in (c) and (d). 
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FIG. 5. The same as in Fig. 4. The pairing symmetry of 
the superconductor is d-wave. 



a) 










2 -2 



-1 

E/A„ 



FIG. 6. The same as in Fig. 4. The pairing symmetry of 
the superconductor is /-wave. 
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FIG. 7. Normalized tunneling conductance (j{E) as a 
function of E/Aq for zq — (solid line), zq — 2.5 (dotted 
line), for different orientations; (a) and (c) for surface orien- 
tation normal to a-axis, (b) and (d) for surface orientation 
normal to 6-axis. The pairing symmetry of the superconduc- 
tor is dxy-wa.ve in (a) and (b), _y2-wa.ve in (c) and (d). 
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FIG. 8. Tunneling conductance cr(-B) in the presence of 
a magnetic field /ibH/ Aq = 0.2 parallel to the a-axis of the 
crystal, as a function of the energy E/ Ao, for surface orienta- 
tion normal to a-axis. The strength of the barrier is z — 2.5. 
The pairing symmetry of the superconductor is (a) p^-wave 
(py-wave), solid (dotted) line, (b) da;-wave, (d^-wave), solid 
(dotted) line, (c) fx-wa,ve, (/^-wave), solid (dotted) line, (d) 
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